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oo 
o 
o 



F. Lechenault 1 , O. Dauchot 1 , G. Biroli 2 and J. P. Bouchaud 3 

1 CEA Saclay/SPEC, URA2464, L'Orme des Merisiers, 91 191 Gif-sur-Yvette, France 

2 CEA Saclay/SPhT, UMR2306, L'Orme des Merisiers, 91 191 Gif-sur-Yvette, France 

3 Science & Finance, Capital Fund Management, 6-8 Bd Haussmann, 75009 Paris, France 



O 

a 

O 

i cj i 

<N 
> 

m 

i— H 

vo 

o 

r-- 
o 



PACS 
PACS 
PACS 
PACS 



64.70.Pf - 
05.40.Ca- 
45.70.Cc - 
61.43.Fs - 



Abstract. - The dynamical properties of a dense horizontally vibrated bidisperse granular mono- 
layer are experimentally investigated. The quench protocol produces states with a frozen structure 
of the assembly but the remaining degrees of freedom associated with contact dynamics control 
the appearance of macroscopic rigidity We provide decisive experimental evidence that this tran- 
sition is a critical phenomenon, with increasingly collective and heterogeneous rearrangements 
occurring at length scales much smaller than the grains' diameter, presumably reflecting the con- 
tact force network fluctuations. Dynamical correlation time and length scales soar on both sides 
of the transition, as the volume fraction varies over a remarkably tiny range (5 (/)/(/) ~ 10 -3 ). 
We characterize the motion of individual grains, which becomes super-diffusive at the jamming 
transition </>j, signaling long-ranged temporal correlations. Correspondingly, the system exhibits 
long-ranged four-point dynamical correlations in space that obey critical scaling at the transition 
density 



Introduction. — As the volume fraction of hard 
grains is increased beyond a certain point, the system jams 
and is able to support mechanical stresses. It has been 
argued that this rigidity transition is akin to phase transi- 
tions in thermal systems, and that the jamming density is 
a genuine critical point [1]. Despite recent efforts [1-11], 
the mechanisms underlying the jamming transition remain 
an open problem. "Jamming" appears to be associated 
with two different notions, not always well distinguished 
in the literature. On the one hand, the glass/jamming 
transition reflects the freezing of structural degrees of free- 
dom associated to topological neighborhoods and the di- 
vergence of the structural relaxation time r a , as observed 
in glass- forming liquids, colloidal suspensions and gran- 
ular assemblies. The common mechanism for the diver- 
gence of r a appears to be increasingly collective rearrange- 
ments [12-15], supporting the idea that the transition is a 
critical phenomenon - albeit of a new kind. On the other 
hand, the rigidity /jamming transition, is the appearance 
of mechanical rigidity [1,3,4,19] for which there are only 
a few direct experimental evidence of the mechanisms at 



play. Contrary to the glass/jamming transiton, it does not 
require the system to be at equilibrium: an assemblies of 
grains will eventually becomes rigid at high enough den- 
sity regardless of the rate of squeezing (an extremely small 
rate corresponds to the "equilibrium route"). This rigid- 
ity/jamming transition may in principle be different from 
the glass/jamming transition as indeed found in mean- 
field glassy models [16,17], in a previous experiment [18] 
and below. 

In this letter, we investigate the rigidity /jamming tran- 
sition in an assembly of horizontally vibrated bi-disperse 
hard disks. In contrast with previous studies [8,10,12,14], 
we use a quench protocol that produces very dense states, 
far beyond the glass density <\> g relative to our experi- 
ment, i.e. such that the structural relaxation time r a is 
much larger than the experimental time scale. However, 
there is a density range <\> g < <\> < (j) a where the strong 
vibration can still induce micro-rearrangements through 
collective contact slips that lead to a partial exploration 
of the portion of phase space, restricted to a particular 
frozen structure. For <j> > a , on the other hand, the sys- 
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Fig. 1: The packing fraction axis and the rigidity transition 
with respect to the glass transition and the fully arrested state. 
Previous work on dynamical heterogeneities in granular pack- 
ings concerned the liquid region, approaching the glass transi- 
tion [8, 10, 12, 14], while the present work focuses on <j> > <fi g , 
the glass transition packing fraction, which for our experimen- 
tal timescale is less than 0.84. 

tern is completely arrested (see Fig. 1). We study how 
these structurally frozen, but still evolving states acquire 
or lose their rigidity as the packing fraction crosses the 
value (j)j ~ 0.842. We (i) find that the collective con- 
tact rearrangements build, around 0j, non trivial diffusion 
properties at length scales much smaller than the grain di- 
ameter, (ii) demonstrate that the associated microscopic 
motion occurs through persistent, but disordered large- 
scale currents and (iii) provide the first direct experimental 
evidence of a simultaneous divergence of dynamical cor- 
relation length and time scales at the rigidity /jamming 
transition 0j, distinct from those of the glass/jamming 
transition, and presumably associated with the underly- 
ing heterogeneous dynamics of the force network. 

Experimental set-up and protocol. — Our ex- 
perimental system is shown in fig. 2. A 1:1 bidisperse 
monolayer of 8500 brass cylinders of diameters d srna u = 
4 ± 0.01mm and du Q = 5 ± 0.01mm laid out on a horizon- 
tal glass plate mounted on three micro-controlled stands 
for accurate level control. Power injection is provided by 
a 1.5 kW motor with 4.77 N.m nominal torque, bond to 
an eccentric rotor and rod system that oscillates the plate 
horizontally in one direction at a frequency of 10 Hz and 
with a peak-to-peak amplitude of 10 mm. The grains 
are confined within a fixed rectangular metal frame while 
the bottom plate moves, so that energy is homogeneously 
injected in the bulk of the grain assembly through fric- 
tion. The cell has width L « 100 d sma u, and its length 
can be adjusted by a lateral mobile wall controlled by a 
fim accuracy motorized translation stage, which allows us 
to vary the packing fraction of the grains assembly - the 
control parameter for our experiment - by tiny amounts 
(5 (/)/(/) ~ 5 10 -4 ) within an accuracy of 10 -4 . The pres- 
sure exerted on this wall is measured by a force sensor 
inserted between the wall and the stage. The stroboscopic 
motion of a set of 1500 grains in the center of the sample 
is tracked with an accuracy of 2.10 -3 <i by a CCD camera 
triggered through TTL electronics by a reflector mounted 
on the rotor and probed by an infrared sensor. Lengths 
are measured in d sma ii units and time in cycle units. 




Fig. 2: Experimental system: one monolayer of disks is vi- 
brated horizontally, while controlling the volume fraction and 
tracking the grains 

In order to reach a steady state analogous to those re- 
ported in [20], i.e. in which macroscopic observables are 
in one to one correspondence with the control parameter 
in a reversible and history independent manner, and pro- 
vide reproducible measurements, we first perform a series 
of compression steps to reach a highly jammed configura- 
tion with (j> = 0.8457, before we stepwise decompress the 
system. The pressure signal is acquired as a function of 
time (both with and without vibration) and is shown in 
Fig. 3. In the compression phase, the pressure relaxes af- 
ter each step but never reaches a stationary state. In the 
decompression phase however, stationary (no aging was 
observed) and reversible (see compression cycles on Fig. 3) 
pressure signals are obtained. However, the structural re- 
laxation time r a is much longer than our experimental 
time. Actually, the neighbors of a given particles typically 
do not change on our experimental timescales, indicating 
that we are focusing on a regime where <\> > <p g . Although 
the system is not fully equilibrated, it appears to be in 
partial equilibrium within the basin of configurations cor- 
responding to a given structural arrangement of the grains. 
The mean pressure, its residual static part (when vibra- 
tion is switched off) and the kinetic part defined as their 
difference are shown as a function of the packing fraction 
in Fig. 3, in the decompression phase. At high packing 
fractions, static internal stresses are important, and com- 
pletely dominate the mean pressure. The grain assem- 
bly, in this regime, is rigidly locked to the frame. As the 
packing fraction is decreased, a kinetic contribution ap- 
pears as the grains unlock from the side walls and start 
making a distinctive noise. We identify these phenomena 
with the rigidity /jamming transition, and thereby obtain 
e [0.8417,0.8422]. Whereas the precise value of (j>j 
might depend on the details of the protocol, the critical 
features observed at the transition are expected to be uni- 
versal. Nevertheless, we note that the value ~ 0.842 
is remarkably close to — 0.8422 reported in [19], using 
different grains and a different experimental protocol. At 
lower packing fractions, the pressure is purely of kinetic 
origin. The grains slide on the sides and moves along with 
the plate, and the static part is zero. 

Experimental results. — We now analyze the evo- 
lution of the dynamics of the grains, namely the mean 
square displacement, the self-intermediate scattering func- 
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Fig. 3: Left: Normal stress on the side of the cell normalised 
by the total grains weight, as a function of time during a typical 
run. Inset: Packing fraction protocol. Right: Mean (blue), 
mean static (red) and mean kinetic (green) constraints as a 
function of packing fraction. Vertical segments crossing the 
red circles indicate the range of fluctuations. Inset: Density 
cycles between the two values of the packing fraction pointed 
out by the arrows indicate a reversible and stationary pressure 
signal. The data drawn in red corresponds to periods of time 
during which we stop the vibration and measure the residual 
static constraint. 

tion, and the dynamical correlations, as a function of the 
packing fraction across the rigidity /jamming transition. 

Ratting particles, identified as grains performing instan- 
taneous jumps much larger than the average motion, ex- 
hibit a strongly anisotropic dynamics, but appear to be a 
rather marginal population - ranging from 1 to 5 percent 
of the grains respectively in the highest and lowest density 
states. All dynamical quantities are hence computed for 
non-rattling particles. As a result, and even though the 
macroscopic vibration is unidirectional, motion is found to 
be isotropic - a feature already observed in [8]. We com- 
pute the mean square displacement for a density (j) and 
lag r: 

^w^/^Eii^h 2 ^^)) (!) 



where N is the number of imaged grains and || Ar^|| (t, r) 
is the square displacement of grain i between times t and 
t + r in the grains' geometric center of mass frame. The 
time average (-) t is performed over 10,000 cycles taken 
1000 cycles after each step in <£. The evolution of cr^(r) 
with the packing fraction is presented in Fig. 4. Note the 
very small values of cr^(r) at all timescales, which is a 
further indication that the packing indeed remains in a 
given structural arrangement for the range of densities 
and time scales we are considering here. Although very 
small, ct</>(t) is well above experimental error bars and is 
not the result of noise, as the following analysis further de- 
mo ntrates. At low packing fractions <fi < <fij, and at small 
r the mean square displacement displays a sub-diffusive 
behavior. At longer time, a diffusive regime is observed, 
at least up to our experimental time scale. As the packing 
fraction is increased, the typical lag at which this cross- 
over occurs becomes larger and, at first sight, does not 
seem to exhibit any special feature for <\> ~ qbj (corre- 




Fig. 4: Left: cr as a function of r for = 0.8402, 0.8407, 
0.8413, 0.8417, 0.8422, 0.8426, 0.8430, 0.8435, 0.8440, 0.8443, 
0.8446, 0.8450, 0.8452, 0.8454. Increasing values of tends to 
reduce ct</>(t) (except around <j>j). The bold line corresponds to 
the packing fraction 0.8417 closest to <j)j - note that it crosses 
the data corresponding to cj) < cj)j. Inset: Typical trajectories 
of a grain during a 20 minute acquisition at packing fractions 
smaller (blue) and larger (red) than <f>j. Right Local slope 
v — <9 log ct0(t)/<9 log (t) vs. t for the same packing fractions. 



sponding to the bold line in Fig. 4). Above 0j, an interme- 
diate plateau appears before diffusion resumes. However, 
a closer inspection of cr^(r) reveals an intriguing behavior, 
that appears more clearly on the local logarithmic slope 
v = d log cr ^{r) / d\og{r) shown in Fig. 4. When v = |, 
the motion is diffusive, whereas at small times, v < \ 
indicates sub-diffusive behavior. At intermediate pack- 
ing fractions, instead of reaching | from below, v over- 
shoots and reaches values > \ before reverting to - from 
above at long times. Physically, this means that after the 
sub-diffusive regime, the motion becomes super- diffusive 
at intermediate times before eventually decorrelating and 
entering the diffusive regime. At higher packing fractions, 
this intermediate superdiffusion disappears: one only ob- 
serves a crossover between a plateau regime at early times 
and diffusion at "long" times (but still small compared 
to r a : as noted above, we only explore the regime where 
<T0 <C d). 

In order to characterize these different regimes, we de- 
fine three characteristic times: ti(0) as the lag at which 
u(t) first reaches 1/2, corresponding to the start of the 
super-diffusive regime, t s d(4>) when u(r) reaches a max- 
imum v*(cj)) (peak of super-diffusive regime), and td(4>) 
where v(t) has an inflection point, beyond which the sys- 
tem approaches the diffusive regime. These characteristic 
time scales are plotted as a function of the packing frac- 
tion in Fig. 6. Whereas t\ does not exhibit any special 
features across 0j, both t s b and td are strongly peaked 
at (j)j. Since superdiffusion is tantamount to long-time 
correlations in the motion of particles, this result shows 
that dynamical correlations are maximal at 0j, reinforc- 
ing the interpretation of <j)j as a critical packing fraction. 

As a decisive signature of the critical nature of </>j, we 
now show that these long-time correlations are accompa- 
nied by the growth of spatial correlations of the dynamics. 
We follow previous works [21], where the four-point corre- 
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Fig. 5: Left: Q a (r) and Right: X4,<i(t) as a function of prob- 
ing length a and delay time r for <j) — 0.8413 < <j)j in a log-log 
color plot. The bold line represents 0.613 o>(r). 

lation G4(r, r), and its integral over space X4( r ) were P ro_ 
posed to measure the fluctuations of the temporal relax- 
ation (see also [12,14]). We characterize the local dynam- 
ics through the correlation function Q a (r) = (Q* (x, T)) tx : 

Qifrr) = -jf -x) exp ("^1^ (2) 

where a is the length scale over which we probe the mo- 
tion. Qo,(t) is akin to the self-intermediate scattering 
function in glass- forming liquids. An example of Q a (j) 
for a specific value of <j) is shown in the left panel of Fig. 5 
as a function of probing length a and lag r. The bold 
line represents a = Aa^r), where A = 0.613 is a density- 
independent constant. It is found to define the crest where 
Qo( r ) decreases from one - when no significant motion has 
occurred - to zero - when all grains have made a displace- 
ment larger then a. We found Q a {^~) to obey the scaling 
relation Q a {r) = Q^^r)/ a) at all (/>, showing that cr^(r) 
rescales the whole distribution of displacements at lag r, 
hence being the only microscopic length scale relevant to 
individual motion at this time scale. The four-point dy- 
namical correlation is defined as: 

G 4 ») = (Ql(x + f,r)Q t a (x,T)) tx - <Q*(f,T))J a . 

'(3) 

The susceptibility X4,a(^) = / drG^ a {r,r) gives access 
to a typical number of grains that move in a correlated 
way when the system evolves over a given time scale r. 
An example of X4,a(^) fc> r the same specific value of <\> is 
shown in the right panel of Fig. 5 together with the same 
bold line as above. This time it defines the crest where 
X4,a( r ) is locally maximum: for any given a, X4,a( r ) is 
very small both when r — > 0, when no significant mo- 
tion has occurred, and at large r corresponding to ho- 
mogeneous normal diffusion. Hence X4,a partly inherits 
the strong scaling property found above for Q a and reads 
X4,a( r ) = h(f) ( r ) X (°~cf>( r )/ a ) where the amplitude en- 
codes the relevant collective information. The global max- 
imum xt(4") is reached for a density dependent lag t*(</>) 
and a probing length a*(0) = Xa^ (r*) for which the dy- 
namics is maximally heterogeneous. Thus a* it is not fixed 
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Fig. 6: Timescales n, t s d, td and r* as a function of the pack- 
ing fraction. The lines are guides to the eyes obtained by best 
power-law fits. The dashed line indicates the packing fraction 
cj)j corresponding to the jamming transition 

from the beginning, as often the case in the recent litera- 
ture on X4, but is rather imposed by the dynamics itself. 
We found that r* is always very close to t s d (at least in 
the vicinity of 0j), thus indicating that dynamical correla- 
tions are maximal when the dynamics is maximally super- 
diffusive. A typical snapshot of the Qa(x, r*) at <pj and 
the corresponding displacement field are shown in the left 
panel of Fig. 7; it reveals that the dynamical heterogeneity 
is in the present case organized in channel currents mean- 
dering between blobs of blocked particles. These events 
occur on time scales such that r* « t s d- Hence the ob- 
served currents are persistent in time, and they build up 
the anomalous super-diffusive motion observed in our sys- 
tem - and also reported in soft disks simulations [26] - 
even though the amplitude of the corresponding displace- 
ments is in the present case a very small fraction of a 
grains diameter. More strikingly, we have found that the 
four-point correlation function G4, a *(r, r*) obeys critical 
scaling in the vicinity of <pj (see central panel of Fig. 7): 

n *A G4,a*(0,T*) / r \ 

G 4 ,a* (r, r ) cx Q ^- J (4) 

with Q{u) = exp(— u -71 ), n « — \ and where £4 (</>), plotted 
in the right panel of Fig. 7, is the length scale over which 
dynamical correlations develop. The value of exponent a 
is found to be smaller than 0.15 but cannot be determined 
accurately. This scaling form, together with the strong 
increase of both £4(0) and r*(0) over a minute range of 
0, is our strongest evidence that the jamming fraction <pj 
is indeed a critical point, where a static pressure appears 
and long-range dynamical correlations develop. 

Discussion. — As emphasized above, the typical dis- 
placements are very small compared to the grain size: 
0>(td) <C d. Rather, these length scales compare to inter- 
particle spacing and typical slip events. The topological 
neighborhoods do not evolve: when <j>> <j> g all experimen- 
tal time scales are shorter than the characteristic struc- 
tural relaxation time r a . The exploration of phase space 
hence probes different contact configurations rather than 
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Fig. 7: Left: Typical snapshot of the Qa(x,r*) at <j)j and the corresponding displacement field (magnified 5 times). Center: 
Rescaling of ln[G4,a* (r, r*)/G4,a* (0, r*)] as a function of ^Jr 1^ for 8 densities around 0j. Right: £4 as a function of the 
packing fraction. The lines are indicative power- laws ~ A±/|0 — </>j| 1//2 consistent with our data, but not an experimental 
determination of critical exponents. 



different structural arrangements [22]. 

The length scale, or wave-vector, at which the dy- 
namics is probed is varied continuously over more than 
two decades (see [23] for a similar analysis in a model 
of glass- forming liquid). This analysis reveals that the 
global maximum xl(0) °f the f° ur point susceptibility is 
reached for length and time scales a*(<j>) and r*(0) such 
that a*(0) « a0(r*). This implies two important con- 
sequences. First, the collective dynamics occur on length 
and time scales selected by the dynamics itself and are not 
arbitrarily chosen. Second, it was shown for glass forming 
liquids that the dynamical susceptibility X4,a obeys the 
following general inequality: X4,a > {9Q a / 9(f)) 2 ((j) 2 ) c [13]. 
Plugging in the scaling behaviors Q a (^; 0) = Q (o> (^) /a) 
and X4,o(^) = ^0 ( r ) X ( cr 0( r )/ a )? one finds [25] that the 
temporal variations of X4,a> i n particular the presence of 
a global maximum, are in fact completely encoded in the 
anomalous behavior of cr^(r), proving the strong connec- 
tion between anomalous diffusion and dynamical correla- 
tions. 

Altogether, the jamming/rigidity is found to be a crit- 
ical phenomenon characterized by the emergence of giant 
dynamical fluctuations at the contact scale, embodied by 
extended and collective slip events. We speculate that 
these dynamical fluctuations mirror the underlying hetero- 
geneous force network and its dynamics, since the contact 
network is more likely to slip in the regions of concentrated 
constraint gradients, i.e. along force chains. The mecha- 
nism by which diffusion sets in is therefore through collec- 
tive contact slips allowed by force fluctuations. The decor- 
relation time td corresponds in this picture to the time it 
takes for the dynamics to achieve a complete renewal of 
the force network. Thanks to these micro-rearrangements, 
the effective compression modulus of the system under vi- 
bration turns out to be orders of magnitude smaller than 
that of individual grains, which can be considered as in- 
finitely rigid. We believe that friction plays an important 
role here, by allowing the system to sustain an external 



pressure at densities smaller than <j) a . It would be inter- 
esting to elaborate on Edwards' ideas [24] and understand 
quantitatively the shape of the pressure vs. density ob- 
served in Fig. 3 which appears to scale as P ^ ((/> — <pj) 
close to 0j, as also reported in [19]. 

Hence, at variance with what is often assumed, the 
rigidity /jamming fraction is not the point <p g at 
which the structure freezes within the experimental time 
scales. In fact, studies of glassy systems within a sort 
of mean field approximation suggest that the ideal glass 
(or glass/jamming) transition is distinct from the rigid- 
ity/jamming transition [16,17,28] which is not either the 
point (j) a at which dynamics is totally arrested, as found 
by us and also reported in [18]: an homogeneous diffusive 
regime is found after sufficiently long times both below 
and above <j)j (at least up to the longest available times). 
Only after further compression will the characteristic time, 
Ti, diverge, leading to a complete arrest of the system. In 
our system, we found <p a > 0.845, a lower bound for the 
random close packing fraction. When <pj < <p < <fi a , even 
though the system is mechanically rigid, the vibration en- 
ables it to explore different micro-configurations. This ex- 
ploration eventually generates uncorrelated particles dis- 
placements on longer time scales, hence diffusive motion 
at a very small scale. What will happen at time scales 
much longer than the experimental ones? Depending on 
the volume fraction, either the dynamics starts feeling the 
frozen structure and subdiffusion sets in again or the struc- 
ture unfreezes and diffusion is maintained for ever. Pre- 
dicting such very long term evolution amounts to address 
the question of the existence of an ideal Kauzmann-like 
glass/jamming transition and its relation to the rigidity 
transition. Such a conceptual issue is obviously far be- 
yond the scope of the present study. 

Following recent work [3], it is tempting to conjecture 
that the correlated currents observed here are related to 
the extended soft modes that appear when the system loses 
or acquires rigidity at <pj. Under the action of a mechan- 
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ical drive the system should fail along these soft modes. 
This picture is also compatible with the Mode- Coupling 
Theory (and some Kinetically Constrained Models), where 
deep metastable states become marginally stable at the 
transition [29,31], with a diverging dynamical correlation 
length [30,32]. Since in our case r a ^> r*, dynamics 
is confined within a unique structural arrangement and 
hence activated events are absent, as assumed in MCT. 
More work is needed to confirm (or disprove) these ideas 
and rationalize the value of the critical exponents a, n 
and v*((j)j) associated to the spatio-temporal correlations 
of these structures, as well as the critical exponents gov- 
erning the divergence of r* and £| at <frj. The link with 
the critical properties at the glass/jamming transition (as 
studied in [14] for looser systems <\> ~ (j) g ) also needs to be 
clarified. 

In this work, we have provided convincing evidence that 
the rigidity /jamming transition is a critical point, where 
long range correlations develop both in time and space. 
The jamming density, 0j, where mechanical rigidity sets 
in, is found to be distinct from the density (j) a where the 
dynamics is fully arrested; we suspect that both friction 
and drive are responsible for this intermediate regime. 
The loss of rigidity results from the emergence of large 
scale heterogeneous superdiffusive currents. Remarkably 
the individual displacements themselves are much smaller 
than the grain size pointing out the role of the force net- 
work rearrangments. We expect that our findings may be 
significant to understand the dynamics and rheology of 
dense granular media, colloids and other soft glassy ma- 
terials, where correlated clusters of mobile particles have 
been reported on several occasions, in particular in con- 
nection with soft modes. It may also well be that the 
seemingly ubiquitous 'compressed exponential' relaxation 
in soft glassy materials [15] is related to the superdiffusive 
motion unveiled here [33]. 
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